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STATISTICAL  INFERENCE'.  IN  MARKOV-RENEWAL  PROCESSES 


By  Daniel  P.  Heyman 


I.  Introduction 

The  usual  procedure  in  operations  research  studies  is  to  model 
a  system  and  then  analyze  system  behavior  from  the  model.  In 
applying  and  testing  these  models,  one  is  faced  with  the  problem  of 
estimating  the  parameters  and  testing  assumptions. 

The  Markov- Renewal  process  has  been  used  by  Jewell  [  5]  to 
solve  a  special  structure  dynamic  programming  problem,  and  it  can 
be  used  to  describe  a  queueing  process  with  arbitrary  service  time 
distributions.  This  paper  is  addressed  to  the  problem  of  finding 
estimators  for  the  parameters  of  the  process  and  to  making  infer¬ 
ences  about  stochastic  service  systems. 

The  inference  problem  for  renewal  processes  is  identical  to 
that  jf  making  inferences  about  random  variables;  much  work  has 
been  done  to  develop  a  large  sample  theory  for  Markov  process  (see 
Billingsley  [  ).  It  will  be  shown  that  the  combination  of  these  two 

methods  respectively  suffice  to  find  maximum  likelihood  estimators 
for  the  Markov-Renewal  process  and  to  show  that  with  large 
samples,  chi-square  methods  can  be  used  to  test  composite  hy¬ 
potheses. 

To  illustrate  the  use  of  these  methods,  a  procedure  for  esti¬ 
mating  the  parameters  of  the  arrival  and  service  time  distributions 
of  a  queueing  process  with  Poisson  distributed  arrivals  and  general 
service  times  will  be  outlined.  Also,  it  will  be  described  how  to 


test  the  hypotheses  that  interarrival  and  service  time  distributions 
are  independent  of  the  number  in  the  queue. 

II.  The  Markov- Renewal  Process 

A  Markov- Renewal  process  (M.R.P.  )  is  a  generalization  of  a 
Markov  process  in  which  the  transition  time  from  state  i  to  state 
j  is  a  sample  from  a  probability  distribution  that  may  depend  on 
both  i  and  j  .  The  first  investigations  of  processes  of  this  type 
were  made  by  Levy,  Smith,  and  Takacs;  in  references  [  7]  and  [  8]  . 
Pyke  summarizes  the  initial  work  and  presents  original  contribu¬ 
tions.  Since  these  two  articles  provide  a  uniform  set  of  notations, 
definitions,  and  theorems,  they  will  be  used  as  a  basic  source.  An 
extensive  list  of  references  will  be  found  in  [7]. 

The  following  notations,  definitions,  and  theorems  will  be 

used; 

^  1 
J  =  the  state  of  the  system  after  n  transitions, 
n 

X  ^  the  length  of  time  the  process  is  in  state  J  ;  i.e.  , 
n  ^  ^  n 

the  length  of  time  between  the  (n)^  and  (n+1)*'  transi¬ 
tion;  the  distribution  of  X  is  not  degenerate  at  zero. 

n  ® 

J  and  X  are  both  random  variables.  The  process  can 
n  n 

enter  m  states,  where  m  is  either  a  positive  finite 
integer  or  infinity. 


If  the  system  is  in  state  i  heading  towards  state 
the  (n+1)**  transition,  then  J  =  i  . 


,  during 


2 


DEFINITION  Z.l;  Q  =  (Q. .)  i*  the  matrix  of  functions 

m 

satisfying  (i)  Q.  .(t)  =0  for  t  <  0  ,  and  (ii)  E  Q. .  (ac)  =  1  , 

"  j=l 

( 1  <  i  <  m  +  1)  . 

DEFINITION  2.  2:  The  m  vector  A  =  (a.,  a^ . a.,  .  .  .  ) 

i,  c  j 

is  the  vector  of  initial  probabilities;  i.  e.  ,  ~ 

DEFINITION  2.  3:  The  (J,  X)  process  is  a  two-dimensional 
stochastic  process  n>0}  that  satisfies 

(i)  Pr  (Jq  =  k)  =  aj^ 

(ii)  Xq  =  0  a.  s.  ,  and 

(2.1)  Pr{J^=k.  X^<x|Jq.  Jj.  X^.  X2'“"*'^n-1' ^n-1^ 

^  Q.  (x)  for  all  x  c  (  oc)  and 
‘'n-l 

1  <  k  <  m  +  1  . 

From  this  definition,  we  can  interpret  con¬ 

ditional  probability  that  the  next  transition  is  from  state  i  to 
state  j  ,  that  the  transition  will  occur  in  t  or  less  time  units, 

given  that  the  system  is  currently  in  state  i 
A  " 

S  =  r  X.  for  n  >  0,  S  is  the  elapsed  time  for  the  first  n 
n  .  i  —  n 

1=0 

transitions. 

LEMMA  2. 1;  The  two-dimensional  (J,  S)  process  is  a 
Markov  process,  and  the  J-process  is  a  Markov  chain.  Let 
P. .  =  Q. .(  oo)  ,  then 

(2.2)  Pr{J^=k.  <  y|jQ,  Jj,  Sj .  } 


a.  s. 


n-1,  k 
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and 


(2.3) 


Pr{  Jn  =  j  I  Jq.  .  '  }  =  Pi. 


Since  the  Q-functions  need  not  be  probability  distribution 

functions,  we  define  S..  =  p. .  ^Q. which  "normalizes"  the 

Q. .  ,  if  p..  ^  0  ,  to  the  probability  distribution  function  of  the 

tj  ij 

transition  time  from  i  to  j  .  If  p. .  «  0  ,  ®. .  ^  U,  (the  unit 

*^11  11  1 


step  at  0).  Let  4».j(x)  =  ^ 


H.(t)  =  Z  O.  .(t)  is  the  probability  density  function  of  the 
‘  j=i 


■  V 


t  dH.(t)  is  the  mean  time  in  state 


time  spent  in  state  i;  rj- 
i  The  mean  transition  time  between  states  i  and  j  is  given 

p  or 

by  b.  =  \  t  d'^.  (t).  Several  useful  relations  are  available 
ij  Jo  ‘J 

from  the  definitions  already  presented. 

Pr{X^  <  X  I  Jq . )  =  Hj  (x) 

n- 1 

Pr{Xn<  x|jo . y  J^}  =  ♦j  J  (x) 

n-r  n 
n 

rj.  =  Z  b..p., 

1 

DEFINITION  2.4:  The  integer- valued  random  variable 
N(t)  ^  sup  {  n  >  0,  t  >  0;  S^  <  t}  is  the  counting  distribution 
for  the  number  of  transitions  in  time  t  .  The  integer-valued  ran¬ 
dom  variable  Nj(t)  is  the  number  of  transitions  into  state  j 
before  (N(t)  +  1)  transitions. 


m 


It  is  obvious  that  N(t)  =  Z  N,(t) 


J 


DEFINITION  2.  5:  Set  n(t)  =  (Nj(t).  N^It) . N.(t).  .  .  .  } 


(j) 


.  is  denoted  by  F. .  in  [  7] . 

ij  "  ij  ' 
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The  stochastic  process  f  r)(t),  t  >  0)  is  called  the  Markov- 
Renewal  Process  determined  by  (m.  A,  Q). 

This  means  that  a  Markov-Renewal  process  is  a  vector* valued 
process  which  describes  the  probability  of  having  n.  observations 
of  state  k  in  N  changes  of  state,  given  the  number  of  states, 
the  initial  state  probabilities,  and  the  probability  law  governing 
changes  of  state. 


LEMMA  2.2:  If  m  <  oo  ,  then  for  all  initial  states  i 
Pr  {N(t)  <  n  V  t  >  0  )  =  1  . 


This  lemma  is  of  great  importance  because  it  allows  one  to  fix 
an  interval  of  observation  and  then  record  information  about  the 
changes  of  state,  and  it  implies  that  N^(t^  <  or  . 


DEFINITION  2.6;  The  Markov  chain  determined  by  the 
J-process  alone  is  called  the  corresponding  Markov  chain,  C.  M.  C. 

Man y  relationships  between  the  M.  R.  P.  andlheC.M.C.  can  be 
shown,  but  it  is  of  primary  importance  that  the  conditional  probabili¬ 
ties,  p.j,  are  completely  determined  by  the  C.M.C.  If  the  times  be¬ 
tween  transitions  from  i  to  j  are  observed,  they  will  form  an 
embedded  renewal  process  with  "failure"  distributions 
There  are  m^  embedded  renewal  process  in  the  M.  R,  P.  ,  and 
these,  the  C.M.C.  ,  and  the  vector  of  initial  probabilities  are 
sufficient  to  describe  the  hence  the  M.  R.  P. 

III.  Limit  Theorems  from  Pyke  [6], 


Define  a  function  f  on  (J  ,,  J  ,  X  ).  and  for  t  >  0  define 

n-l  n  n  — 

N(t) 

W(t)  =  r  f  (J.  ,,  J  ,  X  )  Here,  f{ . )  is  a  function  on 

n  =  l 


-  5  - 


i 

1 


the  transitions  and  the  transition  times.  By  picking  f(J^  X^,  )  = 

X. when  J  ,  =  i  and  J  =  j  ,  and  zero  otherwise,  f(  ■  ,  .  ,  .  ) 
represents  the  transition  times  from  i  to  j  .  If  f(J^  1' "^n’ ^n^  ~ 

Jn  H . )  rerepresents  the  n*^  transition;  hence,  W(t) 

can  represent  the  C.  M.  C.  or  the  embedded  renewal  processes. 

The  random  variable  {U  0  )  is  f(J  ,,  J,  X  ) 

summed  over  the  states  after  the  s  realization  of  state  j  ,  up 
to  and  including,  the  (s  -f  1)**  occurrence  of  state  j  .  = 

m^  ,  is  the  mean  of  the  first  passage  time  distribution  in 

M.  R.  P.  ,  and  is  the  mean  of  the  first  passage  time  distribution 

in  the  C.  M.  C. 

THEOREM  3,  1  (Strong  Law  of  Large  Numbers).  Under  certain 
regularity  conditions. 


W(t)/t  —  m  /m  (a.  s. ) 

m. 

the  limit  being  zero  if  p..  =  oo  ,  and  — ^  does  not  depend  on  j 

THEOREM  3.2  (Weak  Law  of  Large  Numbers).  If 
satisfies  the  Strong  Law  of  Large  Numbers;  and  if  is  finite, 

then 

W(t)/t 


THEOREM  3.3  (Central  Limit  Theorem).  For 


p. .  <  00  , 

JJ 


W(t)  -  E(W(t)) 

rVar{W(t)})*'^ 


L 


N(O.l)  . 


IV. 


tie4i 


Statistical  Inference  in  Markov  Processes  f  1]. 

Consider  an  s  state  Markov  chain  with  transition  probabili 


when  n  realizations  of  the  process  are  observed,  let 
»J 


-  6  - 


f.j  =  the  relative  frequency  of  transition  from  i  to  j  and  f.  = 
the  number  of  transitions  oot  of  tftate  i.  •  .J5«fii»€  thf  random  vjriable 


C..  =  (f..  -  f.p..)/f.  ,  then 

IJ  IJ  iMj  I 


f. 


E{?.j}=0  and  plim;^=p. 


n  ‘  i 


THEOREM  4.1  In  a  stationary,  ergodic,  s  state  Marko/ 


chain. 


I 


(f..  f.p..)^/f.p..  , 

ij  -  I'^ij  I'^ij  *d.-  1 


j  p. .  >0 


where  d.  is  the  number  of  the  non-zero  p..  . 

i  ^ij 


LEMMA  4.1.  2)  f  ln(f, . /f.n, .)  ^x^,  n.  p.  .  >0  . 

^  ij  '  i-  ij'  *s{b-1)’  *^ij 

i.  j 

^  f . . 

The  likelihood  functions  is  L  =  >  p.V  .  To  maximize  InL 

L  »J 


L  J 


subject  to  the  conditions  that 


I 


p^^  =  1  for  each  i  ,  form  the  Lagrangian 


AdnL,  X.)  =  ^  In  p..  X  ^  X.(l  -  ^  p. .) 

i  j 


f.. 


^/VlnL.  X  )  =  0  =>  ^  -  X  =  0 


g^A(lnL,  \.)  =0  r— >  X.  =  f. 


Therefore,  the  maximum  likelihood  estimator  of  p. .  is 


fi,  .=f.  ./f.  and  InL  =  )  f. .  ln(f . . /f.)  . 

*^ij  ij  I  max  ^  ij  ij  1 


J 
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The  statistical  Analysis  of  Markov  chains  can  be  accomplished 
by  applying  theorem  4.1.  A  test  for  goodness  of  fit  of  p^j  can  be 
made  by  substituting  P^j  *  lemma  4. 1  (which 

follows  from  theorem  4.1)  gives  the  Neyman-Pear son  criterion. 


V.  Maximum  Likelihood  Estimators  for  the  Markov-Renewal  Procei 


The  probability  that  a  transition  from  i  to  j  occurs  in 
(x,  X  +  dx)  is  q^j|>;»dx  =  dQ_(x)  =  p.j^.j(x)dx  .  In  observing  the 
system,  1'..  samples  from  .'x)  will  be  taken.  Denote  the 
aggregate  of  these  sample  values  as  x^^  ^  and  let  6^^  = 


fl. 


B^),  n  <«,  be  the  parameters  of  , 


where 


6^  may  be  different  in  all  .  The  problem  at  hand  is  to  find 

the  maximum  likelihood  estimators  for  the  p  's  and  the  vector 


components  of  all  6.^  . 


1 


The  "likelihood"  function  of  the  observations  is: 

f 


L  =  "IXj  [  P  d..)]‘ij 

1 1  d  ^  ^ij  ij'  :j  i;" 


where  d  is  the  set  (i,  j)  for  which  f_  >  0  Taking  log?k- 


rithms, 


InL  =  ^  {f. .  Inp. .  +  f . .  In4> .  ,  B  ..]]  ,  and  maximizing 

Li  'J  u  »J  u  u 


InL  subject  to  the  restraint  that 


of  the  Lag  rang  ian 


=  1  for  each  i  by  use 


A(lnL,  X.)  =y  {f..lnp..+  f..  In4>..(x  e.)  -t- )x.(l  -  Ep. .) 

I  L  ‘J  'J  ‘J  L,  ^  j  iJ 


The  function  L  is  called  the  "liklihood"  function  because  it  is 
asymptotically  equivalent  to  the  Radon-Nikodym  derivative  of 
the  process. 


-  ti  - 


1 


we  get 


■gfr^lnL-  =  0 


AClf'L,  X.)  =  1  -  =  0 


j^^lnL,  X.  )  =  gy— ■  (fyln+. j(5^J  ' .  9  ^  ^ 

ij ,  k  ij ,  k  ■'  ■" 


L. 

The  first  two  equations  yield  p..  = ,  and  the  third  is  the 

tj  tj 

nnaximum  likelihood  estimator  of  the  function  (X  )  considered 

ij  ij  - 

only  «  a  renewal  density  function  assuming  the  values  X  . 

This  analysis  points  out  some  major  properties  of  Markov- 
Renewal  processes;  that  the  C.M.C.  and  the  imbedded  renewal 
processes  formed  by  the  state-to- state  transitions  are  asymptoti* 
cally  mutually  independent,  that  knowledge  of  each  of  them  separately 
implies  knowledge  of  the  complete  process  and  vice>versa. 

A  simple  example  in  computing  maximum  likelihood  estima¬ 
tors  will  illustrate  the  results  of  this  section.  Assume  that 


-  O' 't 

Pr  {  transition  from  i  to  j<t}=l-e  'J.  Then 


f..  f..  -e..-Ex.. 


where  the  summation  is  taken  over  those  x.^  representing  transi¬ 
tion  times  from  i  to  j  .  The  Lagrangian  is 

AdnL,  X.)  =  Inp..  .  f..  Ine,.  -e,.  ^X..).^X,(l-^p,.) 


Equating  the  partial  derivatives  with  respect  to  p^^  and 


to  zero  and  solving  simultaneously  gives  p.^  : 
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a 

57T: 


(InL.  \.)  -- 

I 


{. 

M 

If 

U 


X  =  0  =>  ft.  -  {  I  X.  . 

U  »J  u/  _  U 


The  preceding  result  shows  that  the  maximum  likelihood 
estimator  for  the  exponential  transition  time  parameter  between 
states  i  and  t  is  the  number  of  i  ~*’J  transitions  divided  by 
the  sum  of  the  transition  times.  This  is  the  maximum  likelihood 
estimate  of  the  exponential  distribution  parameter  when  the  phe¬ 
nomena  generating  the  observations  is  not  specified. 


VI.  Estimation  in  Queues  with  Poisson  Arrivals  and  General 
Service  Distributions 

If  the  customers  queued  before  a  service  facility  are  observed 

at  the  points  of  departure,  and  is  the  numoer  in  the  queue  im- 

mediatedly  after  the  departure  (n  =  0,  1,  Z,  .  .  .  ),  with 

defined  as  the  time  between  .1  and  J  ,  ,  the  queue  has  been 

described  as  a  Markov-Renewal  process.  The  validity  of  the  pre- 

ceeding  statement  is  established  in  the  following  manner:  Let 

be  the  number  of  arrivals  during  the  service  time  of  the  nth 

customer,  then  7  -  .1  ,  +  C  -  1  if  the  system  was  not  empty 

n  n- 1  n 

8  t 

aft''r  the  (n  -  1)  customer  completed  service;  if  J  ,  -  0, 

n-  1 

I  =  f.  Letting  7J(x)  denote  Heaviside's  Unit  Function,  the  two 
n  n 

relationships  are  <.  ombined  as 


.7  J  ,  -  U{J  ,)  *  ?  . 

n  n- I  n- I  n 

Since  Poisson  arrivals  have  the  property  that  depends 

only  on  X  .  then  J  depends  only  on  J  ,  and  the  J-process 
n  n  n  - 1 

IS  a  Markov  chain.  The  implicit  assumption  that  service  times  are 
independent  samples  of  a  positive  random  variable  suffices  to 


-  iO  - 


identify  the  imbedded  renewal  processes.  If  the  process  has  been 
operating  for  a  long  time,  the  stationary  probabilities  are  the  initial 
probabilities,  and  by  the  comment  following  the  definition  in  2.6, 
the  sufficient  conditions  to  define  an  M.  R.  P.  are  met. 

The  service  time  distribution  may  depend  on  i  and  j  ; 
indeed  it  is  the  null  hypothesis  that  this  is  not  the  case,  which  may 
be  of  primary  interest  when  making  statistical  tests  of  the  model. 

The  state  spaces  of  the  C.M.  C.  consist  of  the  nonnegative  integers, 
since  that  is  the  range  of  possible  queue  lengths  (number  in  queue 
does  not  include  the  number  in  service). 

Notice  that  the  proof  of  the  existence  of  the  C.M.  C.  is  not 
invalidated  if  the  Poisson  parameters  is  dependent  on  i  ,  the  ^alue 
of  It  ie  well  known  that  Poisson  ariivals  with  parameter 

have  interarrival  times  that  are  expoentially  distributed  with  par¬ 
ameter  X.  .  Observing  the  interarrival  times  when  the  system  is 

A 

in  state  i  ,  we  have  from  Section  Y  that  X.  =  the  number  of  arri- 

1 

vals  divided  by  the  sum  of  the  interarrival  times.  The  chi-square  or 
Kolmogorov -Smirnov  tests  can  be  used  to  test  the  goodness  of  fit  of 
X.  and  the  null  hypothesis  that  X.  =  X,  =  .  .  .  =  X.  (i.e.  .  that  the 
arrival  rate  doesn't  depend  on  the  number  in  queue).  If  the  hypothe¬ 
sis  is  rejected,  alternate  hypotheses  about  the  dependence  of  the 
arrival  rate  and  i  can  be  tested  with  chi-square  or  Kolmorogov- 
Smirnov  statistics. 

In  the  preceding  section  it  was  shown  how  to  find  estimators 

for  and  p..  .  In  this  case,  ^..(X,  6..)  is  the 

ij  ij  ^ij'  ij' 

service  time  distribution  when  i  >  1  .  It  is  the  convolution  of  the 
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interarrival  time  distribution  and  the  service  time  distribution  when 


i  =  0  ,  and  the  data  for  the  estimators  is  generated  by  looking  at  the 
i-to-j  transitions  at  the  points  of  departure  from  the  service  facili¬ 
ty.  If  W(t)  is  taken  so  it  represents  the  i-to-j  transition  times, 
and  applying  the  Central  Limit  Theorem  (Theorem  3.  3),  the  chi- 
square  test  can  be  used  to  test  hypothesis  about  the  service  times, 
such  as  v.j  =  4>^  .  all  the  4).^  are  equal,  or  4)^^  depends  on  i 
and  j  in  a  specified  manner.  For  the  i  =  0  case,  must  be 

known  so  the  estimators  for  4>. .  and  chi-square  test  "theoretical 

ij 

value"  items  may  be  calculated. 

At  this  point  the  interarrival  time  distributions  and  service 
time  distributions  have  been  estimated  and  tested.  If  the  queue  dis¬ 
cipline  is  the  same  as  assumed,  the  model  should  accurately  predict 
q.  ,  the  probability  that  there  are  i  customers  in  the  queue  at 
steady  state.  Assuming  that  observations  are  taken  on  a  steady 
state  queue,  q^  =  the  time  spent  in  state  i  divided  by  the  observa¬ 
tion  time  T  .  In  terms  of  the  M.  R.  P.  , 

n 

qj  =  Nj(t)t,./T  =  Nj(t)  2  •’ijPij/T. 

j=l 

The  estimate  of  the  p.  is  p.  =  , 

U  U  \ 

m 

q.  =  N.(t)  y  b../T 

j  =  l 

This  result  should  agree  with  the  calculated  q;  If  this  is 
not  so,  further  investigation  should  be  made  as  to  the  nature  of  the 
queue  discipline  and  those  characteristics  of  the  process  which  were 
not  included  in  the  model. 
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A  NOTE  ON  DEFICIT,  EXCESS,  AND  SPREAD  IN  A  MARKOV  RENEWAL  PROCESS 

by 

Wliliam  S.  Jewell 

1.  Introduction 

Consider  a  Markov  Renewal  Process  (MRP)  with  a  finite  number  of 

states,  started  in  state  i  at  time  t  =  0  .  Let  N(t)  be  the 

number  of  transitions  in  (0  ,  t]  ,  ^  t  the  time  at  which  the 

last  previous  transition  occurred,  and  >  t  the  time  at 

which  the  next  following  transition  will  occur.  PJandom  variables  of 

interest  are:  (i)  the  deficit,  T_(t)  =  t  -  S,  ^  ;  (ii)  the  excess, 

D  TH  t ; 

‘  spread,  . 

The  purpose  of  this  note  is  to  show  the  distributions  of  these  r.v.s., 
particularly  in  the  limit  as  t  -*  oo  ,  We  shall  use  the  definitions, 

notations,  and  conventions  of  FYKE  ^  vho  has  previously  found 

the  distributions  of  excess,  shown  in  slightly  different  form  as  (3) 
and  (6)  belcw. 


2.  Distributions  of  Deficit,  Excess,  and  Spread 

We  first  find  the  probability  that  the  last  previous  transition 
was  into  state  J  during  the  interval  [t  -  y  ,  t]  ,  and  the  next 
following  transition  is  into  state  k  during  the  Interval  (t  ,  t  +  x]  . 
This  probability  is 


,  X  ;  t)  equals 

P[Zt  =  J,  t  -  y  ^  ^  t  <  ^  t  x  I  Zq  =  i]  , 

for  i,  J,  k  =  1,  2,...,m  ;  O^y^t  ;  0<x^oo;  and  t  ^  0  . 


By  a  straightforward  use  of  the  definitions  of  a  MRP,  this  may  be 


expressed  as 


X,  t)  .  -  x)  -  Qjk(t)lUg(y  -  t) 

(1) 

The  desired  marginal  distributions  of  deficit,  excess,  and  spread  are 
then  easily  found  to  be: 

Dj^^(y;  t)  =  P[Z^  =  J,  JN(t)+i  *  *^1)  i  y  1  Zq  -  i]  0  ^  y  ^  t 

(2)  =  -  Qj^(t)]UQ(y  -  t) 

v-y 

E<‘)(x;  t)  =  P(Z^  .  J,  .  k,  Tj  i  X  I  2^  .  1)  0  <  ,  i  .. 

(3)  •  .  X)  - 

*Io  ‘V'‘  *  -  -  -  «jk'‘  ■  ‘ 

and 
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t)  = 


P[Z^  = 


J,  J 


N(t)+1 


'tq  ^  W  I  Z_  =  i] 


(M 


=  Qj^(w)[Mij(t)  -  M^j(t  -  w)] 

■  J,  tu'*’  0  <  “  i  t 

w  “W 


r  t 

■  Jq  •  t  <  w  ^  00 


The  lower  limits  on  x  and  w  may  be  extended  to  include  zero,  be¬ 
cause  of  the  convention  that  F  (O)  e  Q  (o)  =  G  (O)  =  M  .(O)  =  0  . 

IJ  IJ  IJ 

3.  A  Theorem  on  the  Limiting  Distributions. 

As  In  [3l,  PP-  125^-55,  it  is  possible  to  use  a  renewal  theorem 

due  to  SMITH  to  examine  the  limiting  forms  of  fl)  -  (4)  as  t  -•  oo  , 

whenever  state  J  is  recurrent  and  b  <  «  .  The  result  when  state 

Jk 

j  is  transient  follcws  from  the  fact  that  lim.  M,  (t)  is  finite. 

tr^  jJ 

Paralleling,  then,  the  proof  of  Theorem  1.1  in  [3]  ,  we  may  show  that; 
THEORY 

(i)  If  state  J  is  transient,  then  lim^  D^f^Cy;  t)  =  lim^  E^f^Cx;  t) 

"  . . .  Xr^oo  JK  tr*cc  JK 

liin.j^_^sj^\w;  t)  =  0  ,  for  all  i  ,  k  and  all  w,  x,  y  ^  0  . 

(ii)  If  state  J  is  recurrent  and  b.,.  <  00  ,  then 

— —  — — —  —  JK 
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(5) 


:6) 


t)  = 


and 


t)  = 


"Jk’<'') 


^’) 


=  G  (oo) 
^0 


[1  -  Fjj^(z)]d2  -  w[l 


where  it  is  understood  that  if  G  is  a  lattice  d.f.  then  both  t 

^  J 

and  V,  X,  or  y  may  only  take  on  values  equal  to  multiples  of  its 


U.  Remarks 

In  many  instances,  the  identity  of  the  last  previous  and  the  next 
following  transitions  are  unknown.  Sunning  over  the  Indices  J  and  k  , 

we  obtain 


(8i 


D^^^(yj 


h‘'(z)  dz 

J 
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I 


I 


(9) 


(x) 


D^^^(x) 


and 


(10) 


Hj(z)  dz  - 


where  the  P  .  are  the  stationary  (time-averaige)  probabilities  of 

the  MRP,  and  H.(.)  =  1  -  H  (.)•  Various  conditional  distributions 
J  J 

can  also  be  obtained. 

In  the  event  that  the  underlying  Markov  chain  has  a  single  re- 
cxirrent  class  which  is  positive,  it  is  well  kncwn  that  P  P^  for 

^  J  V 

every  J  in  that  class,  zero  otherwise.  In  other  words,  fornr.!.las 
(8)  -  (10)  become  independent  of  the  starting  rate,  and  provide  generali¬ 
zations  of  the  distributions  of  deficit,  excess,  ar.i  spread  of  renewal 
theory,  (cf.,  for  e.xample,  [1]). 


5.  Applications 

The  above  results  arose  in  the  investigation  of  a  model  of  human 
performance,  where  it  was  necessary  to  determine  how  random  sampling 
of  the  length  of  Job  in  progress  biased  the  actual  dlstricution  of 
lengths  of  Jobs  of  a  certain  type. 

The  results  elIso  have  application  in  calculating  certain  temical 
Z'ewards  which  might  arise  in  Markov-renewai  progranmirg  . 
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NEAR-PARALLEL  CONSTRAINED  OBJECTIVES  IN  INTEGER  PROGRAMS 


by 

David  W,  Matiila 


Some  algorithms  for  the  solution  of  integer  programming  problems 

« 

(e.g.  [l],  [2])  search  for  the  optimal  Integer  valued  solution  In  a 
way  that  necessitates  examining  all  feasible  extreme  point  solutions 
where  the  value  of  the  objective  exceeds  its  value  at  the  optimal 
integral  solution. 

The  following  example  shows  that  this  search  can  be  a  prohibitively 
costly  method. 

Example ; 

Maximize  z  > 

(1  *  +  (1  +  (jMxg  *  (1  +  +•■•+  (1  ♦  (y)  )*„.!' 

subject  to  x^  ■  0,  1,  lul,...,n, 

2x  +  2x^  +  2x  4-..,+  2x  +x  "n, 

125  n-1  n 

where  n  is  ar  odd  integer,  n  =  2k  +  1  . 

Since  n  is  odd,  x  must  equal  unity  in  any  feeisible  integral 
solution.  Observe  that  the  constraint  is  symnetric  in  the  first  n-1 
variables  and  that  exactly  ■  k  of  them  must  be  unity  in  any 

feasible  integral  solution.  It  is  clear  that  z  will  have  a  maxima], 
value,  z^  ,  over  the  integral  solutions  if  the  k  x^'s  with  leu-gest 
coefficients  in  the  objective  are  picked  as  the  unity  valued  • 
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The  optimal  integral  solution  is  therefore: 


\  “  "2  -  "3  =  \  -1 

\^i  ’  V2  “  •••  “  Vi  •  ° 


i=l 

f'iow  for  any  1  ^  <  J  ^  <  . . .  <  <  n-1  ,  1  ^  Jq  1  n-1  , 

where  ^  i«l,...,k  ,  let 

*1  l«d,...,k 

•^i 

0  1 
X  *  “ 

1  2 

x*^  ■  0  m  ^  J  l-0,l,...,k 

n  1 

To  show  that  this  is  a  basic  feasible  solution  observe  that  It 

satisfies  the  constraint  and  that  if  it  were  the  midpoint  of  two  solutions 

p  ,  q  ,  both  p  and  q  would  have  to  have  unit  coordinates  for 

variables  ,  Jg,  ,  zero  coordinates  for  variables  m  , 

i»0,l,...,k  ,  and  to  satisfy  the  constraint  x,  ■  «  .  But  then 

Jq  « 

paq  and  hence  we  have  an  extreme  point  solution. 

The  value,  z'*  ,  of  the  objective  for  this  solution  satisfies 

00  1  k 

z-'  ^k  ^  |-k  ^  ^  (^)  >  k  ♦  -  z° 

i^  1^ 
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4.. 


Hence  this  solution  is  a  basic  feasible  solution  with  a  value  of  the 

objective  exceeding  its  value  for  the  optimal  integral  solution.  But 
2k 

there  are  k(j^  )  ways  of  picking  the  J^  above,  each  yielding  a 
different  z''  .  I.e,,  with  only  fifteen  vaurlables  (nnl5)  ,  24,024 

such  bewlc  feasible  solutions  would  exist. 

The  essential  difficulty  in  the  above  integer  program  is  the 
condition  that  the  objective  is  "nearly  parallel  constrained";  that 
is,  the  coefficients  of  the  objective  are  approximately  proportional 
to  the  coefficients  of  a  constreiint.  The  example  given  here,  where 
all  the  coefficients  in  the  objective  (except  one)  are  themselves 
about  equal,  was  chosen  only  to  highlight  the  difficulty,  and  its  very 
special  structure  should  not  cause  it  to  be  considered  a  mere  lith¬ 
ological  example. 

Indeed  if  any  subset  of  the  coordinates  of  the  objective,  independent 
of  their  own  relative  size,  are  all  nearly  the  same  constant  multiple 
of  the  corresponding  coordinates  of  a  constreilnt  (or  linear  combination 
of  constraints),  all  basic  feasible  solutions  to  the  program  that  differ 
only  among  these  variables  will  yield  approximately  the  same  value 
of  the  objective.  This  cluster  of  basic  solutions  might  all  have  to 
be  searched  before  arriving  at  the  optimal  Integral  solution.  As 
there  are  in  general  many  constraints  and  many  subsets  of  variables  that 
could  be  nearly  pareJ.lel  constraioeii,  this  phenomenon  should  not  be 
dismissed  as  unlikely  in  a  real  problem. 


-  22  - 


references 


[1]  Szvare,  W.,  The  Mixed  Integer  Linear  Programalng  Problem  When  the 
Integer  Variables  are  Zero  or  One.  School  of  Industrial  Administration, 
Carnegie  Institute  of  Technology,  May  7,  I963. 

[2]  ELmaghraby,  Salah  E. ,  An  Algorithm  for  the  Solution  of  the  "Zero- 
One"  Problem  of  Integer  Linear  Programming,  Depeirtment  of  Industrial 
Administration,  Yale  University. 


-  25  - 


REVIEW: 


PROGRA>flENG  UNDER  NONLINEAR  CONSTRAINTS  BY  UNCONSTRAINED 
MINIMIZATION:  A  PRIMAL- DUAL  METHOD 

By  A,  V.  Flacco  and  G.  P»  McCormick,  Research  Analysis  Corporation 
Technical  Paper  RAC-TP-96,  September,  196% 

This  paper  is  concerned  vith  the  mathematical  programming 
problem:  Minimize  f(x)  subject  to  >  0  (i  =  l,...,m)  . 

The  method  of  solution  follows  a  heuristically  presented  idea  of 
C.  W.  Carroll--The  Created  Response  Surface  Technique  for  Optimizing 
Nonlinear  Restrained  Systems,  Operations  Research  9(1961),  I69-I69-- 
in  which  ein  effort  is  made  to  tremsform  the  original  problem  into  a 
sequence  of  unconstrained  minimization  problems  in  order  that 
existing  techniques  (e.g.,  steepest  descent)  might  be  brought  to 
bear.  Carroll  proposed  the  minimization  of  functions 

m 

P(x,  r)  =  f(x)  +  r  ^  W^/g^(x) 

1=1 

where  r  and  the  are  positive  constants  and  the  sequence  of 

minimizations  arises  from  the  stipulation  that  r  «  r^,  r^,,..  where 
(r,  )  is  a  strictly  decreasing  sequence  converging  to  0  .  (The 
function  P(x  ,  r)  is  the  created  response  surface.)  It  was  con¬ 
jectured  that  for  each  r  ,  an  x,  would  be  determined  such  that 

K  K 

llnij^  Xj^  a  x  ,  where  x  solves  the  original  problem. 


-  24 


The  authors  (Fiacco  and  McCormick)  make  the  proposal  rigorous 
by  stating  hypotheses  clearly  and  proving  the  necessary  theorems.  It 
is  shown  that  the  constants  in  r)  may  all  be  taken 

as  1  .  A  considerable  amount  of  computational  experience  is  included. 

It  is  worthwhile  listing  the  hypotheses  which  are  made  in  the 
course  of  the  development 

1.  R°  =  (x  I  (x)  >  0  ,  i  l,...,m}  is  non-enrp^y. 

2.  f  ,  g^(l  "  l,...j,m)  are  twice  continuously  differentiable. 

5.  For  every  real  n'umber  k  ,  the  set  of  feasible  points 

such  that  f(x)  k  is  bounded. 

U.  f(x)  ^  >  "  *  ^’cr  all  feasible  x  . 

r, 

5.  A  method  exists  for  obtaining  on  element  of  R'  . 

6.  f  Is  a  convex  function. 

7.  Is  a  concave  function  ''l  1^,,,  m)  , 

8.  P(x,  r)  is  strictly  convex  for  all  x  e  R'^  and  all  r  >  0  . 

The  proof  of  convergence  of  the  method  to  an  optlrvl.  solution  Includes 
conditions  6  and  7  which  constitute  the  assumptions  for  convex  pro¬ 
gramming.  The  authors  contend  that  conditions  1-5  are  "iict  very 
restrictive."  Certainly  condition  1  Is  necessary  to  the  method, 
and  is  not  unusual.  Without  condition  4,  there  would  be  no  problem, 

■  Fiacco  showed  [Operations  Research  9  (l96l),  I84-I85]  that  if  R° 
is  non-empty,  but  none  of  its  elements  is  at  hand,  t’.en  one  con  be 
determined  by  use  of  the  method  Itself.  Hence,  condition  5  iS  not 
particularly  troublesome.  Condition  8  can  obvdously  be  met  in  several 
ways,  one  of  which  is  by  the  inclusion  of  the  constraints 

^  0  (i  «  l,...,n)  into  the  problem.  Conditions  2  and  5  are 
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probably  more  restrictive  than  the  others.  Note  that  the  latter 
precludes  any  horizontal  asymptotic  behavior  on  the  part  of  f  ,  and 
hence  the  optimal  solution  will  occur  at  a  finite  point. 

It  appears  that  the  sense  in  which  this  is  a  "primal-dual  method" 
is  that  the  successive  minimizations  of  P(x,  r)  yield  feasible 
solutions  to  both  the  primal  and  dual  problems  and  hence  valuable 
Information  on  convergence. 


Richard  W.  Cottle 


REVIEW  : 


La  M^thode  du  Chemln  Critique  -  Application  aax  Programmes  de  Production 
et  d' Etudes  de  la  M^thode  P.E.R.  T.  et  de  ses  Varlantes  (The  Critical 
Path  Method  -  Application  of  the  PERT  Method  and  its  Variants  to 
Programs  of  Production  and  Pleumlng),  A.  Kaufmann  eu^d  G.  Desbazellle, 
Dunod,  Paris,  196^+,  ItiO  pages,  2U  francs. 

This  monograph  provides  an  introduction  to  the  now-classical  methods 
of  time-only  (PERT,  etc.  )  and  cost-time  (CPM,  etc.  )  scheduling 
of  projects,  in  which  a  network  is  used  to  describe  the  partial  ordering 
between  the  Jobs  or  activities  of  the  project.  The  scope  of  the  book 
is  fairly  indicated  by  its  Table  of  ^'ontents,  and  some  of  the  sub-topics: 


Chapter  1.  Graphs  and  Orderings.  A  Graph,  Strict  Order  Relations 
in  a  connected  graph  without  circuits  [cycles]. 
Decomposition  of  a  graph  into  levels,  Finding  an  ordering. 


Chapter  2.  Establishing  a  Reseairch  or  Production  Program.  Re¬ 
presentation  of  a  program  by  a  graph,  Critical  path. 
Float  intervals  and  Job  margins,  Finding  the  critical 
path  [two  algorithms  and  one  numerical  example]. 

Two  examples  of  application.  Jobs  with  random  duration. 

Chapter  3.  Generalization  of  the  PERT  Method.  [4  la  H,  Eisner, 
Opns.  Res.,  Vol.  10,  No.  1,  I962]. 


Chapter  4.  Optimization  of  the  Cost  Economic  Function.  Decreasing 
the  total  cost  of  a  program.  Acceleration  of  a  program 
at  least  cost  [numerical  example].  Arbitrary  variation 
of  operating  cost  as  a  function  of  the  duration.  Linear 
variation  and  the  parametric  linear  program,  Fulkerson's 
[emd  Kelley's]  algorithm,  Deted-ls  of  the  iterative 
procedure.  Example,  Optimization  of  a  program  when 
the  duration  of  the  Jobs  are  random  variables  [not  really]. 


Appendices.  The  p  Distribution. 

Mean  Value  of  a  Critical  Path  [4  la  D. R.  Fulkerson,  Opns. 
Res.,  Vol.  10.  No.  6,  1962J. 

Activity  Subdivision  and  Realization  of  a  Schedule.  [4  la 
T.L.  Healy,  Opns.  Res.,  Vol.  9.  No.  5,  I96I]. 
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As  In  other  works  by  Kaufmann,  the  main  strengths  of  this  monograph 

are  that  it  is  timely,  the  explanations  are  clear,  and  the  authors  do 

not  hesitate  to  give  detailed  numerical  examples.  Certeiinly  the 

exercises  involved  in  seaj'ching  all  of  the  current  literature  to  extract 

the  pertinent  points  was  no  doubt  very  tedious,  and  the  authors  are  to  be 

commended  on  eliminating  much  of  the  sophomoric  philosophy  and  black 

maigic  present  in  most  popularized  accounts  of  the  subject.  The  reviewer 

is  in  complete  agreement  with  E.  Ventura,  who  states  in  the  preface: 

"I  6im  Inclined  to  think  that  [the  critical  path  mathod's] 
success  is,  in  a  large  part,  due  to  the  fact  that,  in  contrast 
to  the  other  methods  of  C.R.,  it  does  not  require  a  previous 
elaborate  mathematical  training,  and  makes  use  of  concepts  which 
ane  fairly  intuitive. " 

and  with  the  authors’  conclusion  that: 


"The  critical  path  method  in  application  favors  'exactness', 
'cledrvoyance ' ,  and  also  'aud£u:ity';  it  is  a  good  example  of 
the  use  of  mathematics  in  management.  "  . 


IV)  this  reviewer,  it  was  interesting  to  see  how  Anglo-Seixon  terms 
become  transposed  into  French.  To  avoid  the  Academic  conflict  of 
"mlnlmer"  versus  "mlnimaliser",  the  authors  use  [MIN]  throughout;  for 
"labelling",  they  use,  variously,  "etlquetage",  "afflchage",  and 
"marquage".  Students  of  graph  and  network  flow  theory  should  be 
warned  that  "arc",  "arete"  [edge],  "chemin"  [path],  "chalne",  and 
"circuit"  eire  sometimes  used  differently  than  in  the  U.S.  literature. 
"Ordonnancement"  ceui  mean  simply  [eirrangement  in  regular  order,  or 
in  0. R.  contexts  can  mean  "sequencing"  or  "scheduling",  or  both. 
Technical  French  is  changing,  also;  in  addition  to  the  obvious 
coinings  "b^tonnage"  [concreteing],  "habillage”  [making-ready ] , 
"compeictage"  [comp€u:ting  ] ,  it  is  also  apparently  proper  to  refer  to 
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"grutage"  [ "creuielng" ],  Sometimes  the  French  is  better;  how  marvelous 
it  would  be  to  have  a  title  like  "organisateur ",  or  "pleinifeicteur 
The  "float"  or  "slack"  aasociated  with  Jobs  becomes  "marge  libre", 

"marge  totale",  or  "marge  certain",  which  makes  a  happy  distinction 
with  the  slack  in  an  event,  which  is  simply  "inter’/alle  de  flottement" 
(except  in  a  marvelous  typographical  error  on  p,  U2,  where  it  becomes 
"intei^^alle  de  frottement 

There  is  sane  uneveness  in  the  presentation.  The  very  first  sentence 
of  the  book  is,  approximately,: 

"Without  abusing  the  theoretical  aspects  of  finding  ordinal 
relations  in  a  graph,  this  concept  being,  as  we  shall  see, 
none  other  than  a  set  of  elements  between  which  oriented  connections 
exist,  we  think  that  it  is  suitable  to  recall  certain  terms 
from  graph  theoi:*y,a^d  to  show  how  to  use  various  procedures  permitting 
one  to  achieve  ordinal  relations". 

Even  assuming  the  authors  wish  to  avoid  philosophy,  it  seems  a  small 
amount  of  motivation  would  be  welcomed  at  this  point  by  the  non¬ 
technical  reader. 

Surely  not  all  the  graph-theoretic  concepts  in  Chapter  I  and  the  set- 
theoretic  notation  in  Sections  2,  10,  l8,  27  ,  28  and  Appendix  II  ao*e 
needed  to  state  the  algorithms,  since  no  proofs  ire  given,  and  the 
examples  presented  do  not  require  it. 

Also,  giving  a  separate  Chapter  to  Eisner's  Generalisation,  and  Section 
20  to  The  Entropy  Function  seen  maladroit,  when  the  sp«u:e  could  better 
be  used  to  explore  cost -time  methods,  explain  parametric  studies, 
discuss  resource  loading,  give  a  list  of  available  computer  codes,  or 
discuss  other  importeint  omitted  topics. 


This  reviewer  also  feels  there  are  sane  pedagogical  liinitations  In 
the  monograph.  For  exaoQxLe,  the  reader  Is  assumed  to  know  about  11' .ear 
prograzaning,  duals,  and  network-flow  theory,  but  this  knowledge  is 
never  used,  since  nothing  is  proven.  In  fact,  the  two  examples  of 
cost-time  solutions  only  show  cases  where  the  Jobs  eure  contiuedly 
compressed  (or  remsiin  fixed)  as  the  project  time  decreases.  Pity  the 
poor  reader  who  tries  em  example  of  his  own,  and  finds  some  Job  being 
extended  (bounded  flow  variable  re-entering  the  basis).' 

As  a  personal  bias,  the  reviewer  feels  that  the  Beta  distribution 
is,  as  usual,  over-onphasized.  The  authors  are  never  quite  clear  as 
to  what  a  reuidom  Job  duration  means  (a  one-shot  "personal  probability"?), 
whether  the  scheduling  is  to  be  performed  before  or  alter  the  Job 
durations  eu:e  known,  what  corrective  action  will  be  taken  if  the 
scheduling  is  incorrect,  and  so  on.  The  eumlysis  of  Section  52  , 
which  purports  to  carry  out  cost-time  optimization  with  random  Job 
durations,  is  in  facr  misleading,  unless  one  resuis  very  carefully 
between  the  lines  to  find  out  the  real  sissumptions  of  the  analysis. 

Thus,  although  the  monograph  provides  a  good  survey  and  clarification 
of  the  existing  literature,  the  feeling  persists  that  the  authors 
have  not  completely  integrated  the  material  with  their  own  experience, 
or  added  research  results  of  their  own. 

Finally,  the  reviewer  had  hoped  that  a  little  chauvanlsm  would  be 
shown  in  discussing  Fench  research  in  the  critical  path  area.  Although 
little  known  outside  France,  the  contributions  of  B.  Boy  and  others, 
beginning  from  graph  theory  and  yet  motivated  by  real  scheduling 
problems,  provided  a  i)arallel  and  coincident  stream  of  independent 
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scheduling  methods,  quite  simllcLr  to  the  American  developments 
described  in  this  book.  (Cf.  ,  for  example,  B.  Roy,  Contribution  de 
la  Theorie  des  Graphes  a  1' Etude  de  Problemes  d'Ordonnacement, 

Proceedings  Second  International  Conference  on  Operational  Research, 
Aix-en-Provence,  France,  September,  I960,  IXinod,  Paris). 

Summary  of  Tasting:  A  faithful  pressing  by  a  well -known  chateau  of 
a  famous  stock  imported  from  California.  Body:  adequate.  Color: 
clear,  some  sediment.  Bouquet:  heady.  Alcoholic  content:  high,  some 
residual  fermentation.  Price:  popular.  Conclusion:  needs  more  maturing. 


William  S.  Jewell 


REVIEW: 


Les  Problfemes  d'Ordonnancement  -  Applications  et  M^thodes  (Scheduling 
Problems  -  Applications  and  Methods),  Number  2  in  Monograph! es  de 
Recherche  Operationelle.  A.F,  I.R.O.,  "by  a  group  of  specialists  led  by 
B.  Hoy",  IXinod,  Paris,  1964,  15I  pages,  no  price  given. 


This  monograph  may  conveniently  be  regarded  as  a  supplement  to  La 
Methode  du  Chemln  Critique  (see  Review,  this  issue);  it  provides 
further  elaboration  of  critical  path  mathod.s,  showing  how  the  theory 
inay  be  extended  somewhat  to  Include  other  problems  of  scheduling,  and 
providing  several  clear,  interesting  examples  of  French  application 
of  -chese  methods.  To  a  certain  extent,  this  monograph  cuiswers  the 
reviewer's  objection  that  La  Methode  did  not  sufficiently  stress 
early  French  contributions  to  the  area  of  project  scheduling.  Part  of 
this  seeming  relucteuice  is,  of  course,  due  to  the  fact  that  our 
emphasis  on  "publish  or  perish"  has  yet  to  reach  the  Continent. 


The  monograph's  r€d.son  d'etre  is  given  in  the  Forward: 


"It  is  observed,  when  scanning  speciall^ed  Journals  of  all  types, 
that  scheduling  problems  are  beginning  to  take  on  increasing 
Importance.  Not  only  the  mathematicians  work  increasingly  each 
day  on  theoretical  plans,  but  even  industrialists  hesitate  less  eind 
less  in  confiding  to  operational  chercheurs  the  difficulties 
they  meet  with  in  their  planning. 

These  works  remaining,  nevertheless,  very  dispersed,  the 
French  Society  of  Operations  Reseeu*ch  [So.F.R.O.]  Judged  it 
desirable  to  assemble  into  a  single  volume,  a  certedn  number  of 
texts,  for  the  most  part  already  published  on  the  occasion  of 
seminars,  colloquia,  or  congresses,  which  would  form  a  sample  of 
French  work  achieved  during  recent  years,  in  the  domain  of 
industrial  applications,  as  well  as  in  theoretical  reseeirch. 

Thus,  we  hope  to  offer  to  the  reader,  in  a  relatively  limited 
number  of  pages,  the  meeins  to  get  an  idea  of  the  current  state 
of  the  entire  subject." 
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Because  of  the  variety  of  the  various  contributions  to  the  monograph, 
we  shall  discuss  eax:h  Chapter  separately. 

Chapter  I.  Fhysignomy  and  Features  of  Scheduling  Problems  (B.  Roy) 

This  introductory  article  describes  the  general  problem  of  scheduling 
"something"  by  decomposing  it  into  "tasks",  Euid  then  choosing  a 
starting  time  and  duration  for  these  tasks.  This  choice  is  limited  by 
three  principle  types  of  constraints,  which  Hoy  calls  "i»tential" 
(PEHT-type),  "disjunctive"  (mutuaiLly  exclusive  use  of  facilities), 
and  "cumulative"  (resource-loading).  The  article  concludes  with  a 
comparison  of  the  two  types  of  graphs  which  may  be  used  to  describe 
the  partial  orderings  of  the  taisks:  "Task-Potential"  graphs  (those  used 
by  Roy  and  co-workers),  and  "Stage  Event-Potential"  graphs,  (the 
familiar  "arrow"  diagram  used  in  the  U. S.A.). 

Chapter  2.  Operations  Research  in  Construction  (P.  I^aud). 

This  Chapter  is  primarily  an  extended  exan^e  of  the  use  of  critical 
path  methods  in  constructing  an  industrial  building  with  five  bays. 

Some  attention  is  devoted  to  the  possible  savings  in  decco^sing  the 
various  "critical"  operations  for  the  entire  building  into  Jobs  for 
the  individual  bays,  thus  reducing  total  time.  !niere  is  also 
discussion  of  Indirect  and  direct  costs  in  building  as  a  function  of 
the  total  Job  duration. 

Chapter  5.  Scheduling  of  Operations  and  Personnel  in  the  "Coffering- 
Tunnels”  (J.  deRosinskl). 

An  industrial  construction  exsaqile  is  described  which  consists  of 
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repeated  use  of  "coffering-tunnels "  in  building  up  successive  stories 
of  the  main  framework  of  an  apartment  building.  The  primary  limitation 
in  using  critical  path  methods  was  the  need  for  "craneing"  during 
"rotation",  when  the  "tunnels"  are  put  in  place  on  the  new  "cell"  to  be 
concreted.  This  Imposes  a  "disjunctive"  constraint,  which  was  solved  by 
heuristic  means. 

Another  limitation  on  possible  schedules  was  the  need  to  maintain  a 
proper  "rhythm"  of  "work  cadence",  i.e.  to  schedule  the  manpower 
shifts,  number  of  men  in  a  given  area,  lunchtimes,  etc.  in  a  realistic, 
repetetive  manner. 

Chapter  U.  Maintenance  Scheduling  in  a  Refinery  (D.  Carr^). 

A  general  discussion  of  the  critical -path  method  as  used  to  schedule 
planned  maintenance  in  a  refinery,  attempting  a  realistic  evaluation  of 
advantages  euid  dlsadveuitages  of  the  method. 

Chapter  5.  Scheduling  (Xitfitting  Operations  at  Chantiers  del '  Atlanticiue 
(Fh.  Guitaxt) 

The  principal  engineer  of  Chantiers  d'Atlantlque  describes  the  problems 
of  scheduling  outfitting  on  the  liner  "Fiance",  carried  out  with  the 
aid  of  the  Soci^t^  d'Ecoixiraie  et  de  Math^matiques  Appliqu^es  (S.E.M.A.  ) 
and  Compagnie  des  Machines  Bull.  In  this  application  the  technological 
sequence  problem  was  secondary  to  the  problem  of  maintaining  a  smooth 
manpower-loading  curve  of  the  various  manpower  "corporations": 

"metalwork",  woodwork",  ventilation",  "piping",  "painting",  "electricity", 
and  "sail -making".  Because  of  the  many  Internal  and  external  constraints 
on  the  "starting  date"  and  "elasticity"  [rate  of  work]  of  each  Job,  the 
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schedule  was  found  by  means  of  a  "segmented"  form  of  linear  programclng, 
which  permitted  rapid  visualization  of  the  manpower  loading  curves, 
and  selection  of  a  desirable  schedule.  Even  though  it  wets  at  first 
thought  possible  to  find  a  "best"  solution  directly 

"....indeed,  it  was  impossible  to  define  a  numerical  indicator 
which  permits  one  to  satisfactorily  classify  the  various  possible 
schedulings  relative  to  each  other.  For,  whatever  the  definition 
of  this  indicator,  there  will  always  correspond  a  well-defined 
value  for  a  given  schedule,  but,  unfortunately,  the  converse  is 
not  true. " 

Chapter  6.  Contribution  of  the  "Oieory  of  Graphs  to  the  Study  of 
Scheduling  Problems  ( B.  Roy ) . 

This  article  is  essentially  the  same  paper  presented  to  the  Second 
International  Conference  on  Operational  Research  at  Alx-en-Provence 
in  i960.  It  is  of  interest,  primarily  because  It  indicates  the  early 
French  developments  in  scheduling  theory,  which  were  suprlslngly 
peumllel,  and  in  some  sense,  better  than  the  early  American  P.E.R.T. 
developments. 

Chapter  7.  Using  the  Fortran-Complled-List-Processlng-Language  in  the 
Solution  of  a  Scheduling  Problem  (P.  Dax*naut  and  G.  Sandler) 
This  cu'tlcle  describes  the  use  of  a  list-processing  language  in 
developing  a  computer  program  for  solving  If!RT-type  problesw;  it  is 
stated  that  this  type  of  language  Is  advantageous  In  retaining  a 
compact  description  of  the  network  In  rapid-access  memory.  The  program 
written  for  an  IBM  709O  handles  ^,000  Jobs,  and  gives  the  starting 
dates  "wedged  to  the  right  and  to  the  left",  le.  both  "earliest"  and 
"latest". 
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Chapter  8.  Scheduling  Problems  with  Disjunctive  CoriStralnts  (Ph.  Tuan 
Nghiem). 

The  last  Chapter  studies  the  problem  of  "dis junctlve "  constraints  in 
scheduling  problems  when  Jobs  belonging  to  certain  sets  are  required 
to  be  carried  out  during  disjoint  intervals  of  ti;  e,  as  for  example, 

Jobs  using  the  same  resource. 

"The  set  of  minimal  schedulings  is  studied  and  a  procedure  is  found 
which  permits  one  to  move  over  the  interior  of  this  set,  which 
has  the  structure  of  an  "arborecence "  [forest], by  passing  from 
one  schedule  to  all  of  its  successors  by  a  simple  change  in  the 
order  of  carrying  out  the  tasks  suojcct  to  vhe  constructive 
constraints. " 


All  in  all,  this  collection  represents  an  interesting  group  of 
articles  written  at  consistently  intelligent,  ye>.  uncomplicated,  level-- 
something  which  is  not  time  of  most  American  articles  on  scheduling.  If 
the  future  finds  the  French  publishing  more  of  their  applied  research, 
we  sha3.1  have  and  veicome  a  high-level  competition  ih  the  field  of  O.R, 


It  is  interesting  to  note  that  this  monograph  resulted  from  a  study 

group  organized  by  the  Socidt^  Francalse  de  Recherche  Op^rationnelle 

(So.F.  R.  0.  ),  but  that  this  group  merged  in  March  196U  with  the  Association 

Francalse  de  Calcul  et  de  Traitement  de  1 ' Information  (A.F.  C. A.L. T. I.  ), 

so  that  the  monograph  has  the  imprimatur  of  Association  Francalse  d* 

t.cr 

Informatique  et  de  Recherche  Op^rationnelle  (A.F. I.h.O. ),  (7  ,  Rue 

ieme  . 

de  la  Chaise,  Paris,  J  ,  France).  Is  this  also  a  trend  of  the  future? 

.quminfl.rv  of  Tasting:  An  interesting  variety  of  Frence  "amuse-gueulc".  very 
suitable  for  a  "a^gustatlon"  with  La  M^thoae  du  Chemin  Critique. 


Willieun  S.  Jewell 


REVIEW; 


Scientific  Decision  Matting  in  Business;  Beadlogs  in  Operation  for 
Nomnathenaticians .  A.  Shuchman,  Rinehart  and  Winston,  1963* 


The  author  describes  his  book  as  "a  collection  of  the  best  writing 
I  have  been  able  to  find  that  uescribes  the  alms,  methods  and  tools  of 
management  science  or  operations  research  without  recourse  to  technical 
jargon  or  complex  mathematical  symbolism, "  and  he  supplements  the  papers  with 
connective  discussion  so  as  to  present  a  unified  whole. 

Whether  the  articles  represent  the  best  selection  of  nontechnical 
papers  available,  the  writer  is  not  able  to  Judge  since  he  has  not  made 
a  study  of  the  literature  with  Shuchnan's  edm  In  mind.  As  a  whole,  the 
reEudngs  do  present  a  cross  section  of  viewpoint,  and  the  book  is  the 
first  comprehensive,  nontechnical  atteiiQ>t  at  examining  the  role  of  oper¬ 
ations  research  In  business,  which  also  presents  a  synthesis  of  many  points 
of  view. 

Turning  now  to  the  specific  articles  in  Shuchman's  book,  Sbuchman 
Is  to  be  coB^lmented  on  the  general  content  and  organization  of  the 
material.  There  are  four  major  parts  to  the  book:  Itot  1,  What  Is  Oper¬ 
ations  Reseeurch;  Part  2,  The  Methodology  of  Operations  Research:  Models 
a -a  Model  Building;  Fart  3,  The  Methodology  of  Operations  Research: 
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Techniques;  and  Part  4,  Some  Applications  of  Operations  Research. 

In  Part  1,  the  first  paper,  by  Dean  R.  Wooldridge,  on  the  scientists' 
Invasion  of  the  business  world  is  not  a  very  penetrating  discussion 
of  this  subject.  It  introduces  the  Operations  Research  Specialist 
as  a  scientist  called  in  to  make  an  obvious  observation  on  a  militeu'y 
operating  problem,  and  gives  the  reader  the  impression  that  operations 
research  is  another  form  of  traditional  management  consulting  -  now  a 
respectable  occupation  for  scientists  who  bring  the  virtues  of  objectivity, 
quantltatlveness  and  a  capacity  for  learning  new  fields.  The  interest 
of  scientifically  trained  people  in  the  operating  and  planning  problems 
of  business  is,  however,  real  eind  slgnlficeint. 

The  next  paper,  by  C.  C.  Herrmann  and  John  F,  Magee,  under  the  title 
of  Operations  Research  For  Management,  gives  the  usual  pco-ty  line  of 
Operations  Reseeirch  as  an  introduction  of  the  scientific  method  into 
business  decisions,  defining  models  and  measures  of  effectiveness  to  assist 
in  choosing  among  alternative  courses  of  action.  Here  Operations  Research 
is  described  as  overlapping  but  not  identical  with  Statistics,  Accounting, 
Marketing  Research  and  Industrial  Engineering,  with  the  latter  field 
dlstinqulshed  from  Operations  Reseeu'ch  not  in  problem  content  but  because 
"industrial  engineers  usually  apply  established  methodologies  to  their 
problems-"  a  somewhat  myopic  observation.  Equally  near-sighted  ere  the 
comments  that  "their  work  is  generally  restricted  in  scope  to  manufac¬ 
turing  activities"  and  "Industrial  engineering  is  not  commonly  character¬ 
ized  by  the  mental  disjipline  and  techniques  of  analysis  that  eire  commonly 


associated  with  the  physical  scientist." 

The  next  three  papers  written  by  Melvin  L.  Hornl  provide  the  best 
and  most  complete  answer  to  the  question,  "What  Is  Operations  Research" 
by  discussing  the  developments  that  make  Operations  Research  possible, 
the  needs  and  opportunities  for  Operations  Research  and  the  basic 
processes  of  Operations  Research.  As  a  trilogy,  this  discourse  argues 
that  business  declslorjs  are  rational,  that  the  process  Involved  Is 
oeccmlng  too  conqplex  for  Judgement  alone,  with  Increasing  Interactions 
between  Its  parts  through  automation,  requiring  more  Integration  quantl* 
tatlvely  In  plaimlng  euid  control,  and  that  business  Is  a  flow  process 
requiring  Integrated  pleuuilng  In  time,  all  of  which  provide  challenge 
and  opportunity  for  Operations  Research  to  provide  appropriate  quanti¬ 
tative  models  for  analysis  and  synthesis  In  decision  making.  Biml's 
discussion  of  the  beislc  processes  of  Operations  Research  follows  the 
usual  pattern  of  defining  the  problem,  model  building,  evaluating 
alternatives  and  Inqpllmentatlon.  It  fails  to  make  clear  the  fasillles  of 
Operations  Research  models  which  can  be  used  In  the  more  Isqortant 
decision  making  areas.  Ihe  last  paper  In  Part  1,  by  Akoff,  Is  a  didactic 
repetition  of  the  other  papers. 

The  papers  of  I^urt  2  on  Models  and  Model  Binding  are  generally 
clear,  useful  and  would  be  good  reading  for  scientifically  trained  people 
as  well  as  students  of  business.  One  paper  by  Robert  8.  Weinberg  on 
The  Uses  amd  Limitations  of  Mathematical  Models  attes^ts  to  give  a  de¬ 
tailed  description  of  the  kinds  of  mathematical  models  used  in  Operations 
Research,  partlculaurly  those  related  to  marketing,  this  discussion 
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faills  to  provide  a  cleELr  definition  of  behavioral  equations  as  probabalistic 
models,  and,  in  fact,  is  misleading  as  to  the  character  of  a  stochastic 
model,  Weinberg's  behavioral  equations  are  simple  statistical  regressions 
without  modeling  of  the  underlying  process,  being  at  most  statements 
about  averages  in  which  there  is  considerable  ambiguity  in  identification 
of  parameters,  and,  he  falls  to  explain  the  hazau’ds  involved  in  extra¬ 
polation  and  mailpulation  known  to  most  econometricians  today,  resort!  ng 
finally  to  Justification  by  stating  that  such  models  "help  us  to  system¬ 
atically  organize  our  data  for  analysis."  T}ie  principal  deficiencies  of 
this  part  of  Shuchraan's  collection  are:  (a)  it  does  not  provide  a 
good  description  of  the  activity  analysis  model  of  linear  programming, 
despite  the  long  paper  by  Vazsonyi  o:.  *he  Uses  of  Matnematics  in 
Production  and  Inventory  Control,  (t)  no  goou  examples  of  stochastic 
models  au-e  given,  such  as  those  used  in  queueing  and  inventory  theory. 

The  last  paper  by  Ackoff  on  Prototype  Models  provides  a  reasonably 
complete  classification  of  Operations  Research  models,  but  is  weak  in 
clarifying  the  quantitative  structures  involved. 

P^t  5  0^^  Operations  Research  Techniques  is  by  far  the  most  complete, 
comprising  half  the  book  in  number  of  pages  and  carryirg  the  reader 
through  most  of  the  recognized  techniques.  Shuchmar.  classifies  techniques 
into  those  for  coping  with  Complexity,  Variability,  and  Lack  of  Informa¬ 
tion,  which  is  not  bad,  adthough  one  wonders  why  he  includes  Factor 
Analysis  under  Complexity  rather  than  Variability  where  he  covers  sta¬ 
tistical  techniques.  Dyriamlc  programming  is  classified  under  Complexity 
rather  than  Variability,  which  does  not  emphasize  the  real  significance 
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of  this  method,  out  one  must  admit  that  aJ.i  eiassification  systems 
create  arbitreiry  compartments. 

Under  Tools  For  Copi.n^  With  Complexity,  there  ame  papers  or:  Math¬ 
ematical  Programming,  Dj'namic  Programming,  Symbolic  Logic  and  Factor 
Arradysis. 

The  leading  paper  on  Mathematical  Programmir.g  by  Henderson  and 
Schlaifer  is  a  good  account  of  why  and  when  linear  programming  may  be 
used  in  a  variety  of  single  period,  static  allocation  situations, 
involving  many  Interactiiv'  vau'iables  not  easily  handled  by  trial  and 

error  methods.  It  does  not  indicate  the  full  potential  of  the  technique 

foi'  multistage  decisions  over-  time,  nor  the  limitations  of  uncertainty 

r-elated  thereto,  and  may  be  regarded  by  the  "Executive'  as  a  computer 

solution  procedure  for  rather  trite  problems  in  the  maze  of  his  decision 
making.  It  would  be  useliil,  also,  to  add  here  a  good  nontechnical 
discussion  of  Legwork  Flow  Ajialysis  extending  into  Critical  Path  Methods 
such  as  PERT,  The  second  paper  on  programmin;'  is  an  anti -climax  and 
does  not  add  much  to  the  discussion. 

The  three  papers  on  d^Tiamlc  programming  emphasize  the  multistage 
character  of  the  tech.nique,  but  treat  only  deterministic  problems  and 
the  discerning  reader  may  rightfully  surmise  that  mathematical  pro¬ 
gramming  could  equally  well  be  used  to  handle  the  class  of  prooxems 
discussed.  Thus,  the  real  significance  .1  P. 

lost,  ruuTiely  to  provide  optimal  policies  under  situations  of  uncertainty, 
having  the  character  of  feed-back  control,  l.e.,  a  policy  stating  do 
this  or  that  depending  upon  the  realized  .alues  of  the  random  variables 
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Involved  In  the  actual  operation  of  the  decision  process  studied. 
Shuchman's  description  of  policy  as  terminology  misses  the  point.  The 
combinatorial  limitations  of  Dyr*anic  Programmir.g  could  have  been  clarified 
in  contrasting  this  technique  with  Mathematical  Programming. 

The  single  pa  er  on  Symbolic  Logic  will  undoubtedly  lose  the 
nonmathematlcal  reader,  although  he  may  grasp  the  significance  of 
propositional  calculus. 

The  two  papers  on  F'  ctor  Analysis  are  overly  simplified  discussions 
of  Analysis  of  Variance  which  may  be  understoo.  better  as  topics  under 
Tools  For  Coping  With  Variability,  and  the  statistical  implications 
need  clarification. 

The  section  headed  Tools  For  Coping  With  Variability  starts  with 
papers  explaining  the  concepts  of  mathematical  probability.  The  one 
by  Warren  Weaver  is  the  best,  conceptually,.  .  Some  confusion  is  intro¬ 
duced  by  the  paper  of  Kurnow,  Glasser  and  Ottman,  which  defines  math¬ 
ematical  probability  as  a  relative  frequency  of  an  event  in  an  infinite 
sequence  of  trials,  rather  than  treating  it  as  an  ideeilized  model  for 
comparison  with  experimental  trials  -  and  the  term  probability  measure  is 
used  without  defining  the  space  on  which  the  measure  is  defined.  These 
matters  of  rigor  are  not  trivial,  because  a  good  deal  of  confusion 
about  probability  arises  from  the  fact  that  the  raeasvire  and  event  space 
are  no^  ''''"efuliy  defined. 

The  next  sub-section  under  Variability  is  devoted  to  Q'ueueing 
Theory  and  consists  of  a  single  paper  written  by  Shuchman  which  gives 
ein  adequate  discussion  of  elementary  concepts,  primai'ily  in  terms  of 
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optimizing/  the  n'onber  of  channels  in  a  service  facility, 
this,  the  two  papers  presented  or  the  subject  of  Decision  Theory  do 
not  attempt  any  discussion  of  sequential  decisions.  Instead  they  are 
cotifined  to  subjective  probabilities  and  criteria  for  optimizing  deci¬ 
sions,  with  little  or  no  development  of  model  structures  in  risk  analyst 
The  reader  will  leave  tnis  section  with  nc  clear  conception  of  a  theory 
of  decision.  Perhaps  the  subject  itself  is  vague. 

The  final  sub-section  of  Tools  For  Coping  With  Variability  contains 
two  papers  on  Gabie  Theory,  defined  in  the  first  paper  by  Shublk  as  a 
study  of  decision  making  in  situations  of  conflict.  Shublk '  s  paper  is 
a  good  discussion  of  games,  strategies  and  the  potential  usefulness 
of  game  theory  in  operations  research. 

The  last  section  of  Part  y  on  Techniques  contains  several  papers 
under  the  heading  Tools  For  Coping  With  Lack  of  Information,  discussing 
mainly  topics  in  Statistical  Inference  and  Simulation  suitable  for  a 
completely  uninitiated  reader. 

Part  4  (Some  Applications  of  Operations  Research)  has  three  classes 
of  papers  describing  operations  research  in  Production  Management, 
Marketing  Management  and  Financial  Management,  characterized  by  Shuchman 
as  "not  fully  representative  of  the  Operations  Research  work  they  are 
intended  to  represent, "  since  limited  mastery  of  techniques  presumed 
for  the  reader  has  prompted  him  to  select  papers  which  he  can  at  most 
expect  to  Impart  "a  deeper  appreciation  and  understanluig  of  the  oper¬ 
ations  research  approach  to  business  problems. "  The  papers  dealing  with 
Production  Management  are  confined  to  economic  lot  sizes  and  Inventory 
control.  Little  or  no  examples  are  given  showing  how  operations 


research  haa  been  applied  to  scheduling,  prx)ductlon  planning, 
distribution,  congestion  in  work  flow,  allocation  of 
facilities,  etc.  Those  on  Marketing  Mamagement  are  equally  well  limited 
in  scope,  discussing  statistical  analyses  for  ineMxirement  of  relationships 
between  sales  and  promotional  effort  to  dollar  volume,  with  one  paper 
on  a  large  simulation  study  to  determine  number  and  location  of  ware¬ 
houses  for  a  national  distributor.  No  studies  of  brand  switching 
or  other  marketing  phenomena  are  given  in  which  stochastic  models  are 
used,  nie  papers  on  Financial  Management  deal  with  elementary  saiqpllng 
and  statistical  estimation  in  accounting  emd  rudimentary  use  of  game 
theory  in  capital  budgeting  and  study  of  mergers  and  acg;uisitlons. 

Ab  Shuchman  has  warned,  his  illustrations  of  applications  of  operations 
research  give  the  Impression  that  uses  of  operations  research  in  business 
eure  at  a  low  level  of  rationalizing  decision  making. 

Where  do  we  stand,  now,  after  this  detailed  review  of  Scientific 
Decision  Medclng  ^  Business?  Are  there  major  eureas  of  operations  research 
omitted?  Evidently  Bellablllty  is  not  regarded  as  part  of  businesr 
decision  making.  Di  some  industries,  liotably  Electronics,  reliability 
of  performance  is  an  over-riding  consideration.  nature,  role  and 

problem  areas  of  oon^ter  systems  is  more  or  less  ignored,  yet  business 
is  deeply  engaged  in  such  systems  for  data  processing,  record  keeping, 
perfoznance  of  clerical  tasks  and  analysis.  Here  Lifomation  Ibeory 
deserves  cotisideratlon.  A.  so  a  more  adeqiuate  discussion  of  control 
processes  in  business  would  be  desirable.  Beplaceswat  and  Maintenance 
of  Equipment  is  another  inqiortant  area  not  receiving  adequate  discussion. 


Strangely,  except  for  elementary  notions  of  costs,  modern  economic - 
theoretic  models  and  concepts  in  operations  research  are  not  seriously 
considered.  Specifically  the  role  of  operations  research  in  evaluation 
of  opportunity  costs  and  planning  capital  investments  deserves  more 
discussion  -  paxticulaj^y  linear  economic  models  ought  to  be  of  interest. 
The  whole  area  of  forecasting,  which  is  of  great  imix)rtance  to  business, 
is  susceptible  of  operations  reseeirch  study  and  certainly  more  useful 
than  abstruse  discussions  of  Decision  Theory. 

It  is  difficult  to  escape  the  feeling  that  Shuchman's  collection 
of  best  nontechnical  writings  is  at  too  low  a  level,  neither  impressing 
the  liard  headed  executive  with  the  accomplishments  of  operations  research, 
nor  stimuli  '  he  brighter  students  of  business.  Yet  the  book  may  be 
the  best  comu.-ndi'un  available, 

Shuchmai:  appeals  to  business  executives  and  students  of  business 
administration  for  serious  attention  to  the  "new  management  technology" 
as  being  "corcerned  with  the  very  core  of  a  manager's  Job  -  decision 
making."  One  may  conjecture  from  the  first  paper  "OperrMons  Research: 
The  Scientists  Invasion  of  the  Business  World"  that  Sh;ichnu  n  correctly 
assesses  that  scientists  are  now  regarding  business  management  as  a 
fertile  field  of  endeavor,  but  he  does  not  pursue  the  full  implications 
of  this  development.  Rather  he  argues  for  the  need  of  managers  and 
students  of  business  "to  make  an  effort  to  understeuid  the  discipline" 
of  operations  research,  which  "need  not  be  that  of  an  "xp)ert,  but  it 
must  be  sufficient  to  enable  the  executive  to  know  when  and  how  the 
expert  can  assist  him.  " 
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Shuchman's  concluding  statement  of  his  introduction: 

"Accordingly,  the  aavent  of  Operations  Research  may  very  well 
mear.  that  managerial  decision  making  is  to  be,  in  the  future, 
less  of  an  art  and  more  of  a  science.  If  tlds  is  the  possible 
significance  of  Operations  Research,  then  students  of  business, 
both  in  school  and  on  the  firing  line,  must  examine  it  care¬ 
fully.  " 

has  implications  which  should  be  elaborated  further. 

In  the  reviewer's  opinion,  managerial  decision  making  and  operating- 
control  of  our  modern  industries  and  -.'r'vices  la  in  the  process  of  passing 
from  professionally-trained  or  experience-educated  students  of  business 
to  operations-rese8Lrch-educa:ed  applied  sclent,  s.  This  is  a  nat^n 
result  of  the  growth  of  ou  implex  modern  teclnology-,  where  a  mere 
famllieLrity  with  operations  research  is  not  sufficient  to  meet  the 
challenge. 

Heretofore,  the  problems  of  business  have  attracted  individuals 
with  a  set  of  talents  peculiar  to  the  quid  pro  quo  activities  of  the 
business  eu-ena:  heuristic,  intuitive  reasoning  aptitudes,  coupled  with 
an  extroverted  personality,  strongly  developed  in  an  ability  to  understand 
or  domiriate  other  humans.  Students  with  these  qualities  have  prospered 
in  the  past,  but  now  face  a  challenge  to  which  they  are  fundamentally 
ill  adapted.  The  technology  of  our  business  society  is  rapidly  being 
changed  into  highly  integrated,  automated  systems,  driven  by  the  machine 
and  not  the  man,  and  requiring  an  enormous  amount  of  data  handling. 
Quantitative  models  developed  by  operatlons-research-sclentists  are  beirg 


integrated  with  data  processing  systems,  so  as  to  euitomaticaLly  perform 
the  large  bulk  of  detailed  decisions  heretofore  regarded  as  being  the 
responsibility  of  management.  This  impersonal  technology  will  perform 
better  than  middle  management,  both  in  energy,  capacity  and  quality  of 
decisions,  and  will  ultimately  obsolete  the  need  and  major  employment 
opportunity  for  business  students,  as  well  as  replacing  a  laj*ge  number 
of  clerical  employees. 

In  time,  executives  may  very  well  be  recruited  from  scientifically 
educated  individuals  trained  in  computer  technology  and  operations  reseeirch, 
with  some  education  in  the  soft  sclmces,  since  it  is  this  background 
viilch  will  be  required  to  understand  the  decision  making  system. 

At  this  time,  it  appears  that  new  curricula  in  industrial  engineering 
will  be  providing  most  of  these  students,  both  because  operations  research 
has  become  a  ma^or  component  of  the  Industrial  Engineering  curriculum 
and  because  t’  e  students  traditionally  study  economics,  business  admin¬ 
istration  and  psychology  against  a  common  background  of  engineering 
technology. 

Thus,  it  seems  that  students  of  business  need  to  re-assess  their 
motivations,  aptitudes  and  educational  goals.  Those  with  a  flair  for 
science  should  seek  deeper  education  in  the  mathematical  sciences  and 
operations  research,  obtain  a  real  understanding  of  computer  technology 
and  de-emphasize  the  many  courses  which  they  normally  take  in  instit¬ 
utional  subjects,  moving  so  to  speak  in  the  direction  of  Engineering. 

The  others  should  realize  that  a  revised  management  role  is  needed  for 
Business  Administration.  Clearly  the  technology  of  Scientific  Management 
falls  short  of  defining  objectives  initiating  new  enterprises  emd 
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organizing  hianan  effort.  Leadership  Is  evidently  needed  for  this 
guiding  function.  Can  students  of  business  undertake  this  responsibility? 
Do  they  need  a  more  liberal  education  rather  than  a  somevhat  competence 
In  Operations  Research?  Are  not  studies  of  behavioral  science.  Industrial 
econcay,  social  objectives  and  public  Institutions  more  Important. 


Ronald  W.  Shephard 


Following  are  some  abstracts  of  dissertations  and  theses  submitted 


to  the  University  of  California  at  Berkeley  on  subjects  pertaining  to 


Operations  Reseairch. 


A  Simulation  Model  of  Airport  Landing  Capacity 


By 

Ernest  Farnsworth  Bisbee  II 
M  S.  ,  IndustrieLL  Engineering,  1961 

ABSTRACT 

The  undesirable  prospect  of  increasing  congestion  on  the  nation's 
airways  has  caused  attention  to  be  focused  more  and  more  on  teraiinal 
areas.  Two  distinct  efforts  to  improve  traffic  flow  have  been  studied 
elsewhere.  One,  by  means  of  careful  control  of  aircraft  separation 
along  flight  paths  leading  to  the  terminal  runways  attempts  to  reduce 
flow  variability,  while  the  other  attempts  to  reduce  service  times 
on  the  runways  by  the  use  of  high  speed  exit  taxlways.  Within  this 
paper,  a  Monte  Carlo  simulation  is  proposed  in  order  to  discover 
relations  between  error  of  control  and  traffic  flow  for  optimally 
located  exits  on  a  runway  used  exclusively  for  landing.  The  simulation 
process  Is  discussed  at  length  and  an  expression  is  found  for  minimum 
sample  size  with  stated  confidence  level. 


Velocity  on  a  Two  Lane  Road 


oy 

Lrnect  F ^rnswort^  Bis  bee  II 
Ph.  D.  ,  Frif-ineering,  19t>2 

ABSTRACT 


The  iimitir^j;  average  velocity  is  found  for  an  Infinitely  fast  test  vehicle 
on  a  two-iane  road  in  terms  of  the  average  vehicle  density  in  each  lane  and 
the  relative  velocity  oetween  them.  The  problem  is  solved  under  the  condi¬ 
tions  of  an  infinitely  long,  straight,  road  without  intersections  on  which 
the  traffic  has  a  fixed  velocity  and  an  arbitrary  distribution  of  si>aclngs 
between  successive  vehicles  for  each  lane.  The  solution  of  this  problem 
leads  to  the  expression  of  a  bound  on  the  fundamental  relation  between  flow 
rate,  density,  and  velocity. 

The  methods  used  are  principally  those  of  renewal  theory  and  the  central 
problem  is  shown  to  oe  a  two-dimensional  general  renewal  process  which  is 
solved  by  a  reduction  to  one  dimension  through  probability  arguments.  Ilic 
form  of  the  solution  is  a  simple  expression  which  Indicated  the  Influence 
of  passing  opportunities  on  the  limiting  velocity.  The  model  emphasizes 
moderate  to  high  densities  and  relative  velocities  between  the  two  streasis. 
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Four  Investigations  in  Operations  Reseeirch 


William  0.  aiattner 
M. 5. ,  Industrial  Engineering,  196I 

i)  Solving  Operating  Problems  with  Electronic  Computers 

ii)  A  Possible  Solution  to  a  Special  Case  of  the  Machine  Loading  Problem 
iii  )Simulation  of  Soedcing  Pit  Operations  at  Geneva  Works 

iv)  Sensitivity  Analysis  Illustrated 
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Scheduling  of  Near  Economic  Lots 


by 

James  Douglas  Coming 
M.S. ,  Industrial  Engineering,  1962 

ABSTRACT 

Economic  lot  sizes  cannot  readily  be  scheduled  as  production  lots  idien 
tvo  or  more  Items  share  the  same  production  equlpaent.  A  aimpls  pro* 
cedure  Is  developed  for  allocatli^  the  svsilshle  aaehixie  time  amoi^  the 
various  Items  so  that  a  feasible  manufacturing  cycle  jresults  vith  pro¬ 
duction  lots  approaching  econoad.c  lot  sizes. 
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A  Delivery-La^  Inventory  Model  with  Qnergency 


By 

Klaus  Hermann  Daniel 
n.^D.  ir.  statistics,  1)61 


ABSTRACT 


In  a  recent  article  E.W.  Barankin  investigated  a  one-period,  one- 
commodity  inventory  model  with  a  one-period-la^  delivery  of  "regular" 
orders  and  with  a  built-in  emergency  provision  which  allowed  the 
possibility  of  placing  an  "emergency"  order  with  instaintaneous  delivery. 
The  quantity  to  be  purchased  in  the  case  of  emergency  was  assumed  to  be 
a  fixed  amount  m  .  The  problem  consisted  in  finding  an  optimal 
regular  ordering  policy  and  an  optimal  emergency  policy  so  as  to 
minimize  a  payoff  function  specified  by  the  structure  of  the  model. 

It  turned  out  that,  if  the  stock  level  lies  below  a  certain  critical 
number,  an  emergency  situation  occurs,  which  causes  an  emergency  order 
of  size  m  .  This  raises  the  initial  stock  level  x  ,  which  can  be 
any  real  number,  to  the  starting  stock  level  x  m  .  If  then  the 
starting  stock  level  is  smaller  than  another  critical  number,  an 
optimal  regular  ordering  policy  calls  for  an  additional  purchase  of 
items. 

The  model  presented  here  differs  slightly  from  the  one  prx>poBed  by 
Barankin.  We  do  not  restrict  the  emergency  quantities  [m^^],  i^,2,...,n 

^  Arrow,  JCarlln  Aud  Scarf,  Studies  la  tbe  Mithrnintl  rn1  Theory  of 
Inventory  and  Production, "  Stanford  University  Press,  I95fl, 
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to  be  constajit,  but  only  require  that  thei-e  is  an  upper  bound  m 
on  the  ra^'s  .  Subject  to  this  restriction  we  search  for  an  optimal 
emergency  policy  of  an  n-stage  problem.  This  modified  approach 
results  in  en  optimal  emergency  pol  icy  described  by  certain  cri  tical 
numbers.  The  critical  number  obtained  for  the  one-stage  problem  differs 
now  from  the  one  establised  in  [2].  The  optimal  emergency  quantity 
to  be  ordered  becomes  a  piecewise  linear  function,  defined  on  the 
set  of  initial  stock  levels. 

Both  models  reflect  certain  assumptions  which  will  be  met  in  many 
practical  problems.  The  model  considered  by  llarankin  reflects  the 
situation  that  an  order  for  immediate  delivery  interferes  with 
scheduled  deliveries  of  items  carried  out  by  a  supplier.  Although 
a  higher  unit  cost  is  charged  for  such  an  order,  compared  to  the  unit 
cost  for  regular  delivery,  the  supplier  may  not  be  albe  or  willing 
to  deliver  instantaneously  any  amount  of  items.  If  the  size  of  the 
order  is  too  small  the  delivery  costs  for  the  supplier  may  offset 
his  revenue  resulting  from  such  an  operation.  If  the  ordered  amount 
is  too  large  it  may  exceed  the  quantity  available  at  the  supplier's 
storehouse.  A  compromise  between  these  two  limitations  is  then 
expressed  by  the  assumption  that  only  a  fixed  amount  of  items  can 
be  delivered  without  delay. 

Our  model  does  not  take  into  account  a  positive  lower  bound  on  the 
emergency  quantities  to  be  purchased.  Any  amount,  so  long  as  it  is 
limited  from  above  by,  say,  m^  will  be  available  in  case  of  emergency. 
This  less  restrictive  condition  on  the  class  of  emergency  orders  will 


-  55- 


result  in  an  optimal  expected  loss  function  (for  any  n-stage  problem) 
which  Is  smaller  than  or  equal  to  the  optimal  expected  loss  function 
of  the  model  stated  In  [2].  But  what  Is  even  more  Important,  Is  the 
fact  that  the  extension  to  the  two- stage  problem  of  the  model  treated 
In  [2]  presents  great  difficulties,  since  the  optimal  loss  function 
in  the  on-stage  problem  is  the  minimum  of  two  convex  functions,  which 
is  not  necessarily  a  convex  function.  However,  the  modification 
outlined  above  reduces  the  optimal  loss  function  to  a  convex  function. 
The  convexity  property  achieved  in  this  way  remains  invariant  for 
the  general  n-stage  problem  &b  will  be  shown  in  the  treatment  to 
follow. 

This  opens  the  way  for  the  solution  of  the  n-stage  problem  of  our 
model  after  a  proper  cleiss  of  cost  functions  entering  the  model  hats 
been  chosen. 
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Optimization  of  the  Design  of  Open  Pit  Mining  Systems 


by 

Beverley  Chew  IXier 
D.  Eng. ,  Industrial  Engineering,  1962 

ABSTRACT 


The  optimization  of  equipment  procurement  and  operational  policy  dominantly 
reduces  mining  costs.  Capital  and  operating  expenditures  for  open  pit  mine 
niaterlals  handling  depend  upon  the  geometry  of,  and  schedule  of  ore  removal 
in  the  ore  deposit;  topography;  mill,  storage,  and  transportation  fajlllty 
locations;  equipment  operating  characteristics  and  conditions;  the  fin¬ 
ancial  cai)abllity  of  the  mining  comp€uiy;  and  the  schedule  of  customer 
demand.  This  analytic  integrated  systems  study  of  open  pit  operations 
and  equipment  selection  guarantees  simultaneous  evaluation  of  Important 
variables,  whose  Independent  study  is  Infeasible. 

A  hyi»thetlcal  mine  model,  formed  by  quantization  of  the  mineral  deposit 
into  geographical  zones,  levels  and  ore  types  and  definition  of  the 
related  haulage  road  net,  is  the  basis  for  materials  handling  reduction. 
Stochastic  euialysis  of  open  pit  loading  and  haulage  forms  the  basis 
for  suboptimizing  equipment  type  selection,  road  grade  determination, 
and  optimal  balance  between  loading  and  haulage  units  for  all  locations 
in  the  mine  model  through  the  study  of  a  network  of  independent,  single¬ 
station  cyclic  queues  with  haulage  units  representing  customers  and  the 
loading  equipment,  the  servers.  Finally,  a  smooth  mine  production  sch- 
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edule  (utilizing  equipment  parameters  previously  suboptimized,  is  developed 
oy  linear  progranmlng  to  optimize  equipment  utilization  and  determine 
capital  equipment  requirements  and  operating  policy  to  meet  customer 
demand.  Production  smoothing,  order  of  ore  removal,  production  capacity 
and  schedule  of  custaner  demand  constraints  are  utilized  in  this  large- 
scale  model. 

Simulation  analysis  shows  that  unit  mining  costs  can  be  reduced  by 
25-50'5t  if  the  capacity  and  economic  life  cf  large  scale  equipment  can 
be  utilized.  It  is  found  by  linear  programming  analysis  that  a  50-75^ 
reduction  in  total  equipment  hour  variations  between  adjacent  scheduling 
periods  may  be  accomplished  within  the  set  of  all  feasible  smooth  mine 
production  schedules. 

All  required  input  data  for  this  emalysis  are  available  or  generatable 
economically,  and  all  models  are  economically  optlmlzable. 


Network  Flows,  Graphs,  and  Integer  Programming 


By 

Ellis  Lane  Johnson 
Ph.D. ,  Industrial  Engineering,  I96U 

ABSTRACT 


The  connection  is  developed  between  graphical  methods  and  the  simplex 
method  of  iineeur  programming  for  solution  of  network  flows,  flows 
with  gains,  and  programming  in  an  undirected  graph.  Determining 
dual  veu'iables,  pricing  out,  and  resolving  degeneracy  are  shown  to 
simplify  for  phase  I  of  the  simplex  method  applied  to  the  network 
flow  problem,  euid  the  resulting  eilgorithm  corresponds  closely  to  the 
labeling  procedure  for  network  flows.  Phase  I  of  flows  with  gains 
does  not  simplify  as  much,  but  a  special  case,  programming  in  an 
undirected  graph,  does  simplify  just  as  much  as  in  network  flows,  and, 
hence,  the  primal -dual  method  has  as  much  advantage.  A  solution  of 
8UI  integer  program  is  given  for  phase  I  of  programming  in  an  undirected 
graph.  The  problem  includes  some  previous  work  on  degree-constrained 
subgraphs,  matchings,  and  coverings.  Alternating  paths  play  a  key 
role. 

Matching  stnd  covering  in  a  graph  extends  to  an  important  class  of 
integer  programs:  packing  and  covering  in  a  family  of  sets.  A 
discussion  is  given  of  two  special  cases  of  the  packing  problem: 
coloring  a  graph  and  finding  a  maximum  set  of  independent  vertices 
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of  a  graph.  The  general  packing  problem  Is  reduced  to  the  problem 
of  finding  a  maximum  set  of  Independent  vertlcej.  No  special  device 
for  integer  programming  is  given  fur  this  class  of  problems,  but  a 
refinement  of  Gomory's  cutting  plane  method  for  general  Integer 
programs  Is  given.  The  algorithm  is  based  on  the  simplex  method 
with  multipliers.  Lexicographical  ordering  Is  replaced  by  an  aid 
hoc  perturbation  technique,  and  a  linear  programming  subroutine 
is  used  to  assure  finite  convergence  and  to  speed  convergence. 
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An  Investigation  of  the  Burn  In  and  Related  Problems 


By 

Michael  J.  Lawrence 
M.  S.  ,  Engineering- Industrial ,  I96U 


ABSTRACT 


Two  problems  involving  the  derivation  of  bounds  on  a  distribution 
with  a  decreasing  failure  rate  (a  D. F.R.  distribution)  are  presented. 

1.  Given  that  an  item  has  a  decreasing  failure  rate,  sharp  upper 
and  lower  bounds  on  the  burn  in  time  to  achieve  a  specified  mean 
residual  life  are  presented.  The  bounds  rely  only  on  the  D.  F.  R. 
assumption  and  a  knowledge  of  the  first  moment  emd  a  percentile 

of  the  failure  distribution. 

2.  An  early  estimate  of  the  five  year  survival  proportion  (conmonly 
called  the  five  year  cure  rate)  is  of  great  interest  in  assessing  the 
value  of  a  treatment  for  a  mortal  disease  such  as  cancer.  Assuming 
that  the  distribution  of  time  to  death  is  D.F.  R.  and  assuming  a 
knowledge  of  the  mean  and  a  percentile,  sharp  upper  and  lower  bounds 
on  the  five  year  survival  proportion  are  obtained. 

In  addition  some  bounds  on  the  hazatrd  rate  and  density  of  a  D.  F.  R. 
distribution  ore  stated  sind  proved. 
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Stochastic  Prograoming 


By 

S.  M.  Slnha 

Ri.  D. ,  Statistics,  1962 


ABSTRACT 


This  investigation  is  concerned  vith  a  linear  programming  problem, 
where  the  coefficients  of  the  objective  function,  the  constraint 
inequalities  and  the  resources  are  random  variables.  Such  a  problem 
is  known  as  a  stochastic  programming  problem.  The  linear  programming 
model  for  such  a  case  however  has  no  meaning  and  it  is  necessary  to 
formulate  a  new  model  to  deal  with  such  cases.  The  uncertainty  aspect 
of  the  problem  suggests  that  it  can  only  be  solved  probabilistically 
euid  hence  one  reasonable  formulation  of  the  stochastic  programming 
problem  requires  that  our  activity  levels  should  be  such  that  with 
a  certedn  preassigned  (high)  probability,  the  total  quantities 
required  for  each  item  should  not  exceed  the  available  quantities 
and  at  the  same  time  should  guarantee  a  maximum  objective  with  a 
preassigned  (high)  probability. 

The  usual  methods  described  in  the  literature  for  att€u:klng  this  type 
of  problem  are  not  applicable  because  of  the  non-dlfferentlabllity  of 
expressions  involving  square  roots  of  positive  semi -definite  forms. 
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Project  Planning  By  Decomposition 


Shallendra  C.  Parikh 
M,S.,  Industrial  Engineering,  Feb.  I965 

ABSTRACT 

This  study  considers  "critlcaj.  path"  networks  which  ai'e  used  for  the 
planning  and  scheduling  of  projects  that  consist  of  well  defined 
sequences  of  individual  activities.  When  the  number  of  activities 
is  large,  it  becomes  difficult  to  prepare  a  network  in  the  high¬ 
speed  memory  of  a  digital  computer.  Also,  in  the  cases  where  there 
are  two  or  more  independent  projects,  which  are  weeikly  inter-related 
by  common  activities,  the  problem  of  efficient  scheduling  of  all  the 
projects  becomes  quite  difficult. 

This  thesis  presents  a  method  to  "tear"  or  "decompose"  a  project 
network  into  several  subnetworks,  schedule  the  subnetworks  and  then 
to  put  the  subnetworks  back  together.  Thus,  it  is  possible  to  prepare 
subnetworks  as  individual  units,  to  store  the  subnetworks  in  the 
memory  of  the  computer  and  to  schedule  and  revise  (updace)  the  projects 
as  needed. 

A  computational  ailgorithm  is  first  given  for  time- only  (PERT)  networks 
then,  a  computational  algorithm  is  given  for  a  cost-time  (CPM)  network 
of  project  subnetworks.  The  later  algorithm  is  a  generalization  of  an 
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algorithm  due  to  Fulkerson,  in  order  to  handle  plecevise-llnear, 
convex,  cost-tljne  curves  for  each  euitivity.  Also,  necesscury  conditions 
for  solution  of  continuous,  convex,  coat -time  curves  for  each  activity 
are  presented.  In  the  appendices,  an  efficient  time-only  algorithm 
and  a  node  re-orderlng  algorithm  are  given. 
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Aiial;,“tical  Models  for  Attack  and  Control  of  Wildland  Fires 

George  Merriman  Parks 
f^i.D. ,  Engineering,  1965 


Decisions  concerning  the  dispatch  of  suppression  forces  to  attack  eind 
control  a  wildland  fire  must  be  based  on  the  physical  characteristics  of 
the  fire  and  the  effectiveness  of  the  forces  seta  against  it.  The 
oehavior  of  a  fire  is  affected  b,  the  fuel,  topography,  and  weather 
conditions  in  the  immediate  area,  eind  the  effectiveness  is  determined 
by  the  type  of  force,  the  training,  and  the  equipment.  There  is  a 
dollar  tradeoff  between  maintaining  large  and  costly  suppression  forces 
to  keep  fires  small  and  employing  small  forces  but  allowing  fires  to 
become  large  and  cause  extensive  damages.  These  decisions  have 
traditionally  been  based  principally  on  qualitative  criteria  because 
of  the  lack  of  established  quantitative  relationships  to  describe  the 
behavior  of  free-burning  fires  and  the  effects  of  various  types  of 
suppression  activity  thereon.  The  objective  of  this  study  was  to 
mathematically  describe  these  relationships  to  enable  finding  the 
size  force  such  that  the  sum  of  these  suppression  and  damage  costs  is 
minimum. 

Three  deterministic  analytical  models  were  developed  to  describe  the 
basic  tactics  of  fire  fighting-  direct,  indirect,  and  parallel.  The 
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size  force  necesseiry  to  keep  total  costs  at  a  minimum  was  determined  in 
each  case  under  various  ass'jmptions  concerning  fire  growth  and  suppression 
activity.  This  development  was  extended  to  include  sensitivity  analyses 
of  the  results,  risk  models  where  decisions  must  be  made  under  uncertainty, 
the  allocation  of  limited  forces  to  simultaneous  fires,  and  least-cost 
rei.:forcement  policies.  'I>ie  relationship  of  elapsed  time  between  the 
detection  of  a  fire  and  the  start  of  suppression  activity  and  the  size 
of  force  required  for  control  at  minimum  cost  was  also  investigated 
within  the  framework  of  each  of  these  models.  'T-^o  stochastic  models  were 
suggested  for  determining  optimal  strategies  to  employ  in  certain  fire 
situations. 

Comparison  of  actual  suppression  and  damage  costs  for  the  calendar  year 
1959  on  the  Plumas  National  Forest  with  the  results  predicted  by  one 
of  the  above  deterministic  models  indicated  that  signlficaj't  dolleu-  savings 
can  be  achieved  in  total,  fire-caused  costs  by  increasing  the  size  of 
the  suppression  force  and/or  reducing  the  ti "e  required  for  Initial 
attack.  The  Implications  of  these  results  for  fire  planning  and  budgeting 
were  discussed  and  significant  areas  for  further  research  in  fire 
modelling  were  indicated. 

Various  statistical  analyses  were  performed  on  data  from  approximately 
14,000  actual  wildland  fires  to  determl  le  historical  value'  for  some  of 
the  fire  growth  and  suppression  force  effectiveness  parameters  employed  by 
these  models. 
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Recovery  from  Interrupted  Service:  A  Queueing  Model 


By 
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ABSTRACT 


This  thesis  discusses,  in  physical  and  mathematical  terms,  the  effect 
that  a  temporary  interruption  of  service  has  on  a  servicing  system 
where  the  steady  state  mode  of  operation  is  such  that  queues  are  small. 

The  mathematical  model  employed  assumes  Poisson  arrival  and  exponential 
service  distribution.  Time  dependent  expressions  are  derived  which 
describe  the  behavior  characteristics  of  such  a  system  following 
interruption  of  service.  Specifically,  it  is  possible  to  predict 
a)  the  average  number  of  units  in  the  system,  b)  the  probability 
that  a  late  arrival  will  encounter  delay  greater  thaui  a  specified 
time,  and  c)  the  averag'’  delay  that  a  unit  entering  the  system  will 
experience.  These  predictions  are  possible  for  the  case  in  which  the 
condition  of  the  system  at  termination  of  interruption  is  not  known 
as  well  as  the  conditional  case  in  which  queue  length  following 
interruption  of  service  is  known.  The  ability  to  predict  such 
performance  characteristics  is  of  considerable  value  in  both  the 
initial  design  and  the  routine  operation  of  such  systems. 
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ABSTRACT 


i 


The  solution  of  a  two-stage  linear  program  under  uncertainty  can  be 
obtained  by  solving  an  equivalent  convex  programming  problem.  We 
characterize  the  solution  set  of  this  equivalent  convex  programming 
problem  as  well  eis  the  functional  to  be  optimized.  Since  we  seek  an 
here-and-r.ov  decision,  this  solution  set  and  the  functional  are 
expressed  in  terms  of  the  first  stage  decision  variable.  If  the 
constraints  of  the  initial  problem,  or  the  probability  space  of  the 
random  vao'iable  of  the  two-stage  linear  program  under  uncertainty 
possess  certain  properties,  then  the  solution  set  of  the  equivalent 
convex  program  is  a  convex  polyhedron.  In  this  case,  if  the  problem 
cannot  be  solved  using  existing  solutions  methods  (e,  g.  linear  programming 
with  upper  bounds,  quadratic  programming)  one  ceui  use  an  algorithm 
developed  in  the  last  section  of  this  thesis. 

The  constraints  imposed  on  the  here-and-now  decision  can  be 
divided  in  two  parts:  first,  the  fixed  (or  technology)  constrsd-nts  where 
.10  uncertainty  is  involved;  second,  the  Induced  constrsdnts  which  limit 
our  first  stage  alternative  to  those  which  allow  a  feasible  solution 
to  the  second  stage  whatever  be  the  values  assumed  by  the  random 
elements  of  the  problem.  The  fixed  constrsdnts  are  expressed  sm 
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lineeLT  equations  and  inequalities  in  terms  of  the  here-and-now  decision 
variable.  The  induced  constred-nts  are  usually  not  available  in  that 
form  and  one  of  the  greatest  difficulties  one  encounters  when  solving 
a  linear  program  under  uncertainty  is  to  determine  if  a  given  decision 
satisfies  or  does  not  satisfy  the  induced  constraints. 

When  there  are  no  induced  constraints  one  says  the  problem  is 
"complete";  e.g.  pure  inventory  models,  transportation  problems  with 
uncertain  demand.  This  problem  leads  to  an  equivalent  separable 
convex  programming  problem.  If  the  random  elements  of  the  problem 
assume  some  specific  distributions  functions,  then  the  equivalent 
convex  program  is  reducible  to  programming  problems  for  which  efficient 
algorithms  exist.  If  the  random  elements  are  discretely  distributed, 
then  the  equivalent  convex  program  is  a  linear  program  with  upper 
bounds;  if  the  random  variables  are  uniformly  distributed,  then  it 
is  a  quadratic  program.  In  general  (when  the  random  variables 
are  not  discrete  or  uniformly  distributed)  we  suggest  an  algorithm 
for  which  an  experimental  compute^  program  has  been  written. 


-  69  - 


